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ABSTRACT 

Bolometric corrections are used in quasar studies to quantify total energy output based 
on a measurement of a monochromatic luminosity. First, we enumerate and discuss 
the practical difficulties of determining such corrections, then we present bolometric 
luminosities between 1 jjtm. and 8 keV rest frame and corrections derived from the 
detailed spectral energy distributions of 63 bright quasars of low to moderate redshift 
(z = 0.03 — 1.4). Exploring several mathematical fittings, we provide practical bolo- 
metric corrections of the forms L iso = ( \L\ and log(Lj SO ) = A + B \og{\L\) for A = 
1450, 3000, and 5100 A, where Li SO is the bolometric luminosity calculated under the 
assumption of isotropy. The significant scatter in the 5100 A bolometric correction can 
be reduced by adding a first order correction using the optical slope, a\ :0pt - We rec- 
ommend an adjustment to the bolometric correction to account for viewing angle and 
the anisotropic emission expected from accretion discs. For optical/UV monochro- 
matic luminosities, radio-loud and radio-quiet bolometric corrections are consistent 
within 95% confidence intervals so we do not make separate radio-loud and radio- 
quiet corrections. In addition, we provide several bolometric corrections to the 2-10 
keV X-ray luminosity, which are shown to have very large scatter. Separate radio-loud 
and radio-quiet corrections are warranted by the X-ray data. 
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1 INTRODUCTION 

Bolometric luminosity is the total energy per second ra- 
diated by an object over all wavelengths in all directions. 
Quasars are notable for emitting their substantial radia- 
tion over the entire electromagnetic spectrum. Th e shape 
of a q uasar spectral energy distributions (SED), see I Wilkes] 
( 2004) for a review, depends on the physical properties and 
structure of the quasar. Specifically, these are the proper- 
ties of the central s upermassive black hole, the surrounding 
accretion disc (see iKoratkar fc Blaea 1 19991 for a review of 
emission from this region), the existence of a dusty torus 
and its geo metry, dust heated by starburst regions in the 
host galaxy l|Netzer et al. I [20071 ). and any radio structures 
that may exist. These physical structures differentiate them- 
selves by radiating in different wavelength bands. In the ra- 
dio, quasars look like kiloparsec-scale relativistic jets that 
can extend spatially beyond the host galaxy; in the infrared 
dust of various temperatures from a parsec-scale torus; and 
in the optical/ultraviolet sub-parsec-scale accretion disks. 
As a result, bolometric luminosity and black hole accre- 
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tion rate (expressed as Lboi/L e dd), both fundamental quasar 
properties, are challenging to determine. 

In practice, bolometric luminosity is difficult to measure 
for several reasons. Firstly, quasar SEDs cannot be observed, 
due to a quasar's broad wavelength emission, without data 
from multiple telescopes. This makes constructing a full SED 
expensive in terms of telescope time; a detailed SED can use 
observations from ten or more telescopes, sometimes using 
more than one instrument on a single telescope. 

Despite this issue, quasar SEDs have still been com- 
piled. Usin g observati o ns fro m 12 telescopes and 16 in- 
struments, lElvis et alj (L1994) constructed an atlas of 47 
quasar SEDs and used them to determine bolometric cor- 
rections, which have been used extensively since. There 
have been other recent investigations of qua sar SEDs as 
well (|Shang et all 1201 ll ; iRichards et al.l 120061), those that 
focus on specific wavelength r e gimes ( e.g.. iMullanev et al.l 
201ll; iNemmen fc BrothertorJ I20T0I; iNetzer et al.1 120071 : 



Vasudevan fc Fabian! 120071 : IShang et al.1 120051)." those that 



focus on specific clas ses of objects (e.g., Chen et all l201ll : 



iGallagher et al.ll2007l). a nd also those of individual objects 
(e.g.. IZheng et al.ll200ll ). Despite these contributions, the 
lElvis et al.l (1 19941 ) mean SED and bolometric corrections are 
still often used. 
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Even when data from multiple telescopes are compiled, 
there are gaps in the electromagnetic spectrum that remain 
unobserved. Gaps in the extreme ultraviolet (EUV) and the 
near-infrared (NIR) to a lesser extent, are problematic when 
calculating bolometric luminosity. The emission in these gap 
regio ns has not been observed and can, at best, be modeled 



Shakura fc Simvaev|[l973l ; iMatthews fc Ferlandl [l987j ; 



i Korista et al. || l997 j ) o rl rn^d from sp ectral properties at 
other wavelengths (e.g.. iKruczek et al.ll201lh . Another ap- 
proach is to breach the gap by simply connecting the data on 
eithe r side with a straight line in log-space (e.g-. lElvis et al.l 
Il994h . 

Quasars are variable in many wavelength regions. Be- 
cause emission at different wavelengths originates in differ- 
ent physical structures, different spectral regions vary on dif- 
ferent timescales. Optical emission can vary on a timescale of 
months to years whereas X- ray emission varies on a shorter 
timescale (|Hook et al.lll994l . and references therein). Build- 
ing a consistent, individual SED requires at least quasi- 
simultaneous observations on these timescales in spectral 
regions that suffer from variability. 

Quasar emission can be significantly contaminated by 
light from the host galaxy, particularly in the near and far- 
infrared. The NIR is important for bolometric luminosity 
calculations and, while emission from the quasar usually 
dominates, the NIR sees a drop in quasar emission and a 
peak in emission from the host galaxy. In order to study 
emission intrinsic to the quasar, SEDs must be corrected 
for host galaxy c ontamination. This is usually done (e.g., 
IShang et al.ll201ll ). but does add a source of uncertainty to 
bolometric luminosity calculations. 

Quasar emissio n is anisotropic (e.g., 

iNemmen fc Brothertonl |2010| ). Anisotropy may vary 
from object to object and can depend on wavelength. An 
improved understanding of disc anisotropy can lead to a 
correcti on for average viewing angle to the optical/UV con- 
tinuum. INemmen fc Brothertonl (1201 0T) u s e the theoretical 
accretion disc models of iHubenv et al.l ((2000) to derive 
bolometric corrections and explore their dependencies 
on accretion rate, black hole mass, and viewing angle. 
IHubenv et all (|2000l ) model the accretion disc in the range 
10 14 - 10 17 Hz as a thin disc that includes geometric and 
relativistic effects. Due to the lack of any contribution from 
an X-ray corona, these models do not provide adequate 
emission at X-ray energies as compared to observed SEDs. 
Despite the shortcomings of accretion disc models, they 
are the best guide for the kinds of corrections that may be 
necessary as a result of assuming isotropy. 

These problems have been overcome to varying extents 
to allow calculations of bolometric luminosity to be made. 
In common practice, bolometric luminosity is estimated by 
scaling from a monochromatic luminosity. Diffe rent scaling 
facto rs exist, e.g. at 5100 A: Lj BO = 1 3.2 \L\ (|Elvis et al.l 
1 1994). Li., n = 9AL a <|Kaspi et al.ll2000l ). L i3 = 1.5 x 9 XL\ 
ilShang et all 120051 ). and L lso = 10.33 AL A l|Richards et all 
2006). Note that we differentiate between the bolometric lu- 
minosity calculated under the assumption of isotropy {Li so ) 
and true bolometric luminosity (L(, ;), which likely differ. 

There is a need for new bolometric corrections based 
on up-to-date SEDs with a new methodology. Pointed ob- 
servations (e.g. I Elvis et afll 19941 ) give more detailed individ- 
ual SEDs than those built from survey data l|Richards et al.l 



120061). though s a mple sizes are typically smaller. However, 
the lElvis et ail l| 19941 ) SEDs require u pdating and those 
built from more re cent observations (e.g. iNetzer et al.|[2007l ; 
IShang et al1l2005r ) lack broad wavelength coverage. In addi- 
tion, the common as sumption of i sotrop y is not a good one. 

The 85 SEDs of Shang e t al.l (|201ll ). intended as an up- 
date to those o f lElvis et al.l "ifl994). are based on pointed ob- 
servations from the current generation of telescopes. These 
SEDs combined with a new methodology provide an oppor- 
tunity to handle the above issues in a thoughtful, consistent 
way to calculate the next generation of bolometric correc- 
tions. 

This paper i s organized a s follow s. Section 2 describes 
the data and the IShang et al.l l|201ll) SEDs. Section 3 con- 
tains a description of the bolometric luminosity calculations. 
We derive bolometric corrections in Section 4 and discuss 
them in the context of previous work. In Section 5 we cal- 
culate X-ray bolometric corrections. Section 6 summarizes 
this investigation. 

We adopt a cosmology with Ho = 70 km s _1 Mpc -1 , 
= 0.7, and fi m = 0.3. 



2 DATA AND SAMPLE SELECTION 
2.1 The Full SED Sample 

We derive empirical bolometric corrections from th e N1R- 
to-X-ray continua of 63/85 of the lShang et all (|201ll ) SEDs. 
This atlas has a total of 85 objects from three different sub- 
samples which are described briefly below. 

• The 'PGX' subsample contains 22 of 23 Palo mar- Green 
(PG) quasa rs in the complete sample selected bv lLaor et alj 
1 19941 . Il997l ) to study the soft-X-ray regime. This subsample 
is UV bright and has z ^ 0.4. The optical-UV region is 
covered by UV spectra from Hubble Space Telescope (EST) 
and quasi-simultaneous ground-based optical spectra from 
McDonald Observatory. 



• The 'FUSE-HST' subsample has 24 objects, 17 
of which come from the Far Ultraviolet Spectroscopic 
Explorer (FUSE) AGN program (lKrissll2000l ). This is a het- 
erogeneous, UV-bright sample with z < 0.5. T he SEDs for 
this s ubsample have quasi-simultaneous FUSE (M 00s et al.l 
120001 ). HST, and Kitt Peak National Observatory (KPNO) 
observations. 

• The 'RLQ' subsample includes nearly 50 quasars 
originally assembled to study orientation; all members of 
the sample have similar extended radio luminosity which is 
thought to be isotropic. The SEDs have quasi-si multaneous 
HST and McDonald or KPNO observations. See I Wills et~ai1 
(1995) and lNetzer et al.l l|1995l ) for additional details on this 
subsample. The blazars originally included in this sample 
are excluded in the SED atlas because of their variability 
due to synchrotron emission from a beamed jet. 



X-ray coverage comes from Chandra, XMM, and 
ROSAT observations that are available in the literature. 
When data were available from multiple sources, Chandra 
and XMM are given preference over ROSAT due to their 
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broader energy coverage and higher sensitivity. NIR cover - 
age comes from 2MASS photometry l|Skrutskie et al.l f2006). 
The 2MASS point source catalog has 79 objects in the sam- 
ple. The sample also has FIR and radio coverage, though 
those data are not used in this study. 

We used the NIR through X-ray (rest-frame 1/xm to 
8 keV) region of the SED for calculating bolometric lu- 
minosity. The data in this region of the SED were col- 
lected between 1991 and 2007. AGN variability, which can 
be a concern at these wavelengths, does not appear to be 
a significant issue in these SEDs. Optical-FUV data were 
taken quasi-simultaneously (within weeks) and near-to-mid- 
infrared emissio n, which arises from a size scale on the order 
of parsecs (e.g., iRaban et all [2006) or less for the hottest 
dust, will vary on a longer timescale of months to years. We 
do not see any obvious discontinuities in emission that might 
indicate variability. 

2.2 Corrections 

In order to assure that the calculated bolometric luminos- 
ity of the quasar does not include emission from the host 
galaxy, it is important to remove any contaminating emis- 
sion from the host. The effects of local intervening Galactic 
gas must also be accounted for. I n the process of construct- 
ing the SEDs. lShang et all |201lf ) applied two corrections to 
the data. 

First, at optical-to-FUV wavelengths the SEDs suf- 
fer from Galactic extinction from dust. This is re- 
moved with the empiric al mean extinction law of 



Cardelli. Clayton, fc Mathis] l| 19891 ) using the dust maps of 



Schlcgel et al 



Il998h 



Second. IShang et al.l (|201 lh corrects for host galaxy 
contamination at NIR and optical wavelengths, although 
according to those authors, host galaxy contamination 
to the AGN light is not large for their UV/optically 
bright quasars. The host galaxy contribution is maximized 
compa red to the AGN in th e H-band for redshifts near 
~0.5 (jMcLeod fc Riekd Il995h s o H-band p h otom etry is 
used to determine host fraction. IShang et al.l (|201ll ) make 
observations with IRTF or H ST or collect them fro m 
the literature for 28 ob jects (jMcLeod fc Riekd Il994al lbl: 
iMcLeod fc McLeodl200ll ). 

For objects with IRTF observations, host fraction 
is determined for the H-band using the procedure of 
IMcLeod fc Riekd [|l994ah . In this procedure, a standard star 
is observed directly before or after the target. The stellar 
profile, or point-spread function (PSF), is then scaled to the 
innermost pixel of a one-dimensional surface brightness pro- 
file for the quasar-plus- host. The fraction of the PSF neces- 
sary to make the quasar-minus-PSF profile turn over in the 
center is then subtrac t ed to leave a good estimate of host 
emission. IShang et all |201ll ) make a minimum and maxi- 
mum subtraction of the PSF to estimate an uncertainty in 
this procedure of a few percent, although we note that this 
method seems to over-subtract and is less trustworthy than 
subtractions using HST data. 

For objects with HST observat ions, two-dim ensional fit- 
ting is possible. IShang et alj |201ll ) use galfit (|Peng et al.l 
120021 ') to fit several galaxy models in addition to the PSF. 
The resulting best fit is then used. Uncertainty in this 
method is a few percent or better. 



For objects without H-band observations the 2MASS 
PSF magnitude is taken to be the magnitude of the AGN 
only and the aperture magnitude is considered total magni- 
tude of the AGN plus the host. This method is used for the 
object shown in Fig. [2] 

Once host fraction has been determined in the H-band, 
removing the host contamination is straightfo rward. An 
ellipti cal galaxy template of NGC 0584 from iDale et al.l 
(2007), scaled in the H-band, is then used to determine host 
fraction in the J and K bands and those magnitudes are 
corrected accordingly. Fig. [2] shows the appropriately scaled 
elliptical galaxy template that is used to determine host frac- 
tion for the J and K bands. 

Host galaxy contamination can also affect emission at 
long wa velengths of the opt ical spectra, although to a lesser 
extent. IShang et al l |201ll ) try to remove host contribu- 
tion here when extracting the spectra. They use different 
aperture sizes to ensure that the host contribution is unde- 
tectable in the final spectrum. 



2.3 The Bolometric Luminosity Sample 

We used a subsample of the SEDs for calculating bolomet- 
ric luminosities. Of the total 85 objects in the SED atlas, 
only 63 have the wavelength coverage appropriate for this 
study. The following objects were excluded from the sam- 
ple due to insufficient X-ray coverage: 3C 175, 3C 288.1, 4C 
30.25, 4C 64.15, B2 1351+31, MC2 1146+111, MRK 506, 
PG 1103-006, PG 1259+593, PG 1534+580, PG 2214+139, 
PG 2251+113, PKS 2216-03, and TEX 1156+213. After the 
above objects were removed from the sample, the following 
lacked sufficient NIR coverage and were also excluded: 3C 
110, 3C 186, 3C 207, 4C 12.40, B2 1555+33, B2 1611+34, 
MC2 0042+101, and PKS 0859-14. This leaves 63 objects 
with full wavelength coverage for calculating bolometric cor- 
rections. Table Q] lists sample members. 

The combined sample consists of 23 radio-quiet (RQ) 
and 40 radio-loud (RL) quasars. The RQ quasars are all 
from either the PGX or FUSE-HST samples and are lower 
redshift (z < 0.5), while the RL quasars come primarily from 
the RLQ sample and more than half have higher redshifts 
(z > 0.5). The RL quasars have an average luminosity about 
6 times higher than those that are RQ. Both RL and RQ 
quasars span approximately 2 orders of magnitude in lumi- 
nosity. While not statistically complete or well-matched, the 
combined sample is comprised of subsamples with specific 
selection criteria. The 'PGX' subsample may be complete, 
although it is not representative of the general quasar pop- 
ulation so we advise caution when using corrections derived 
here on more general samples. The high flux limit of the PG 
sample means that t hese bright objects are also uncharac- 
teristically blue and IJester et al.l (|2005h specifically points 
out that the PG sample is biased against redder spectra. 
The 'FUSE-HST' subsample was selected randomly among 
objects known to have high UV fluxes. The 'RLQ' subsam- 
ple was selected to have similar extended radio luminosity, a 
quantity thought to be independent of orientation. In order 
to find large numbers of these objects it can be necessary to 
go to higher luminosity and redshift. 

We compared our sample in observed parameter space 
to the Sloan Digital Sky Survey (SPSS) D ata Release 7 
(DR7) quasar catalog (jSchneider et al.ll2010l ). Fig. [1] shows 
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Figure 1. The region of observed parameter space occupied by 
the 63 objects in this sample and the SDSS DR7 quasar catalog. 
SDSS quasars are small gray points, SDSS R,L quasars are dis- 
tinguished with small red points, RL quasars from this work are 
blue squares, and RQ quasars from this work are open black cir- 
cles. The bolometric corrections derived here may lose accuracy 
for objects that occupy another region of parameter space. 



the observed SDSS g magnitude vs. z. We calculated ap- 
proximate SDSS g magnitude for our objects by converting 
the flux at the center of the g-band to a magnitude with the 
SDSS flux zero point. The radio loudness parameter, log R* , 
is estimated for SDSS objec ts with FIRST dete ctions at 20 
cm using the prescription of IStocke et a l. (1992) (note that 
we substituted g for B). The bolometric corrections derived 
here will be most accurate for objects that lie in the region 
of parameter space occupied by this sample and may be less 
reliable when used on dissimilar objects. 

While the majority of this sample was selected to be 
bright in the UV, there may be intrinsic reddening in some 
objects. It can be difficult to differen tiate between intri nsic 
reddening and an extreme SED fe.g. JZheng et al]|2001r ). In 
light of this, we did not exclude objects at this point. 



3 DETERMINING BOLOMETRIC 
LUMINOSITY 

To know the bolometric luminosity exactly, it is necessary to 
observe a quasar at all wavelengths and from all directions. 
Since we are restricted to one line of sight and can only ob- 
serve parts of the electromagnetic spectrum, we must make 
some assumptions in order to calculate bolometric luminos- 
ity. These assumptions, and our method, are outlined below. 



3.1 Empirical Isotropic Luminosity 

The measured observed frame flux is converted to an ob- 
served luminosity under the assumption, which we will re- 
examine later ( H4.4J1 . that emission is isotropic. 



- J obs,u i ls 



F 



d L sin 8 d9 d(j> 



Lobs,v obB = 47T d L F bs,v 



(1) 

(2) 



where d_t is the luminosity distance, F bs,v oba is the observed 
flux density, and v ba is the observed frequency. The inte- 
gration between Equation [T] and Equation that results in 
the 4-7T coefficient is possible only under the assumption of 
isotropy. Then the observed luminosity is converted to rest 
luminosity via the equation: 



L - 



bs,u obs 
1 + 2 



47rd| F obs> „ o 

1 + 2 



(3) 



where v re3 t is the rest frequency. Bolometric luminosity is 
then calculated by integrating L rest ,i/ rcst over frequency, 
where specific consideration should be given to the limits 
of integration. 

Emission at IR wavelengths is thought to result from 
optic al-UV photons reprocessed by surrounding hot dust 
fe.g.. lBarvainisiri987l ). Under the assumption of isotropy, it 
is important not to include reprocessed photons. Assuming 
isotropy requires that the optical-UV photons that we ob- 
serve are characteristic of radiation from the accretion disc 
in all directions. We count these photons for the first time 
when we integrate the optical-UV over all space. IR photons 
are likely optical-UV photons originally emitted in another 
direction that have been reprocessed and re-emitted along 
our line of sight by a AGN-heated dust. We count these 
photons a second time if we integrate the IR portion of the 
SED. 

To accurately estimate the bolometric luminosity un- 
der the assumption of isotropy, the IR region should be 
excluded from the integration. There has not been agree- 
ment in the field on this point: some (e.g. lElvis et al.lll994l : 
lRichards~e t al. 2006) i ntegra t e the entire SED whereas oth - 
ers (e.g. iNetzer et all 120071 ; iNemmen fc Brothertonl l2010h 
integrate only the region in emission at optical-through- X- 
ray wavelengths (i.e. the Big Blue Bump). Here we integrate 
the SED from 1 ^m to 8 keV (i/ = 3x 10 14 - 1.93 x 10 18 Hz). 
Note that we assume that there is no significant energy con- 
tribution at very high energies (> 8 keV). This assumption 
is made perforce because the data do not extend past this 
limit. The region of integration is shown in Fig. [2] for a typ- 
ical quasar SED. Energy output at radio wavelengths is in- 
significant compared to emission from the accretion disk so 
it can be excluded without significantly altering the inte- 
grated luminosity. 

An alternative to this method is to use the integrated 
mid-infrared (MIR) emission and co vering fraction as a 
measure of bolometric lumino sity (e.g.. ISpinoglio fc MalkarJ 
ll989l : lEdelson fc Malkanlll986h . This method has several ad- 
vantages since the MIR is less susceptible to reddening and 
likely much more isotropically emitted than the Big Blue 
Bump. In particular, such a correction would be valuable 
for reddened and Type 2 objects. It would not be without 
its own issues, however. Corrections for starlight would be 
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necessary at the shortest wavelen gths and also for st arburst- 
heated dust in the host galaxy (jNetzer et alj 120071 ) at the 
longest wavelengths, depending on the limits of integration. 
Additionally, depending on redshift and available facilities, 
the MIR bump may be more difficult to observe than the 
Big Blue Bump and also has gaps in data coverage, although 
they may be easier to fill. Finally, while the MIR emission 
may be more isotropically emitted, the anisotropy of the op- 
tical/UV emission that acts as a source for MIR reprocessing 
must be accounted for. 



3.2 Integrating in the NIR 

In order to integrate the SED, the NIR region around 
Iog(f) ~ 14.25 — 14.75 required special consideration. Data 
coverage here was patchier than the optical/UV (2MASS 
photometry versus a spectrum) and in some objects the 
2MASS data did not reach the 1 /im limit of integration 
in the rest frame. 

In the NIR we used a power law interpolation between 
2MASS data points and then, when necessary, extrapolated 
out from the longest wavelength 2MASS point to provide 
data all the way to the 1 pirn limit of integration. Note that, 
because we chose a hard 1 /im rest frame limit of integration 
but had a range of redshifts, we integrated different data in 
each SED (e.g. in the object shown in Figure [2] we did not 
integrate between the H and K bands or beyond because the 
1 /jm limit had already been reached). 

In order to quantify interpolating between the 2MASS 
points and extrapolating beyond, we replaced these inter- 
polation/extrapolation po wer laws with the F v oc v 1 ^ 3 
l|Shakura fc Sunvaevl ll973T ) accretion disc model extrapo- 
lated out from the optical spectrum to 1 (im and recalculate 
bolometric luminosity. The v 1 ^ 3 model is reached by inte- 
grating a multiple blackbody accretion disc over infinite ra- 
dius and is appropriate for the long wavelength emission of 
the disc (i.e. the NIR). Both the v 1 ^ 3 and power law inter- 
pol ation models are shown in Fig. [2] 

iKishimoto et al.1 l|2005h find good agreement between 
the v 1 ^ 3 model and the observed emission from quasar accre- 
tion discs. Their result is primarily based on PG 1425+267, 
an object in our sample, as well as two other broad line 
quasars of similar redshift to our sample (several tenths). 

We normalized the u 1 ^ 3 model between 5000 and 5300 
A, a region relatively free of emission lines, and which should 
have little contribution from hot dust emission. Emission at 
these wavelengths may be contaminated by residual host 
;alaxy emission, even after the correction of IShang et al.l 
2011 ). We may have introduced some small uncertainty as 
a result, but this is not larger than is introduced by imple- 
menting the v 1 / 3 disc model. 

We found that the difference in bolometric luminosity 
between our power law interpolation and the v 1 ^ 3 disc model 
is small, on average —2 ± 2 per cent. Bolometric luminos- 
ity resulting from v 1 ' 3 model is plotted against bolometric 
luminosity resulting from the power law interpolation in Fig- 
ure [3] 



3.3 Integrating in the X-ray 

The X-ray region of the SED also required special measures 
to be taken before it could be integrated. There were two 



issues in this region: in some objects we did not have data 
out to our 8 keV limit of integration and there was a large 
gap in data coverage for all objects between the FUV and 
X-ray around \og(v) ~ 15.75 - 17.25. 

In some objects we had to extend the X-ray power law 
to higher energies to reach the 8 keV limit of integration. 
Observations in this region were a magnitude or more down 
from the peak of the SED, so changing this limit did not 
significantly change the integrated luminosity. 

Our region of interest in the SEDs contained a large 
gap in data coverage: one between the FUV and X-ray. 
We bridged this gap with a linear interpolation in log vL v 
space. This power law appeares t o be a good des cription 
of this region in some cases (e.g., lLaor et ail 1 19971 ). Fig. 
shows the power law interp olati on as well as models from 
iMatthews fc Ferlandl (|l987i ) and iKorista et all fll997l '). de- 
scribed below, over-plotted for one representative object. 

While a straight-forward power law did seem to connect 
FUV and X-ray data well in some cases, we recalculated 
bolometric luminosity using two other models to help quan- 
tify the level of uncertainty that can result from selecting 
different models in the EUV. 

Both models were normalized to the data near 1300 A, 
a region also relatively free of emission lines. In cases where 
the model did not intersect the X-ray data, the X-ray data 
were extended to longer wavelengths until they intersected 
the model. 

The first model is the SED of IMatthews fc Ferlandl 
l|l987h . which they derive from continuum observations 
when avai lable and otherwise infer f rom the emission-line 
spectrum. IMatthews fc Ferlandl (|l987f ) model the region be- 
tween the FUV and X-ray with several power l aws. 

The second model is that of lKorista etli] l|l997l ). They 
model the UV region with a modified power law of the 
form F„ cc j,-°- 5 e - WfcT^t ^ thermal bump). At 912 A 
and shorter wavelengths, the power law F v oc v~ x is used. 
We used the SED quoted by those authors as typical for an 
average quasar: T cu t = 10 6 K, E pea k — 44 eV. 

The per cent difference in bolometric luminosity caused 
by using the IMatthews fc Ferlandl l|l987h model instead of 
a power law interpolation is 41 ± 21, where the given un- 
certai nty is the standard deviation. Using the lKorista et all 
l|l997h model caused a percent difference in bolometric lu- 
minosity of 33 ± 20. Visual inspe c tion s howed that in many 
cases, the IMatthews fc Ferlandl {l987) and IKorista et al.l 
l|l997h models may not be a good description of the SED; 
in some objects these models over or under-estimated the 
X-ray data by as much as an order of magnitude for over 
an order of magnitude in frequency. In a few objects these 
models always over-estimated the data by approximately 0.1 
dex. 

Bolometric lu m inosit ies re sulting from the 
IMatthews fc Ferlandl l|l987l ) and IKorista et ail (|l997f l 
models are each plotted against bolometric luminosities 
resulting from the power law interpolation in Fig [3] This 
figure illustrates how significantly our ignorance of emission 
in this FUV- X-ray ga p impac ts our b olometric lumi n osities ; 
IMatthews fc Ferlandl l|l987h and IKorista et ail (|l997f ) 
X-ray models tend to result in much larger bolometric 
luminosities than a power law interpolation with lots of 
variation between objects. 

It has been suggested that the intrinsic SED contin- 
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Figure 2. Data in the NIR to X-ray region of the SED for 4C 19.44 shown by solid black line with the region of integration denoted 
by dashed gray verti c al lin es. In the NIR the i/ 1 / 3 disc model is shown as a dark red, dashe d-dotted line. In the FUV-soft- X-ray, the 
Matthews & Fcrland (1987) model is shown as a blue dashed line and the Korista ct al. (1997) is shown as an orange dotted-dashed line. 
The power law interpolation s that we u l timat ely empl oy are shown as so lid red lines. The solid green line shows the elliptical galaxy 
template for NGC 0584 from lDale et al] d2007h used in lShang et al] l|201lf ) for host galaxy subtractions to the 2MASS photometry. The 
host template is scaled to 15%, the host galaxy fraction for this object, of the 2MASS H-band magnitude. Wavelength and energy are 
given on the top axis. 



uum may be different from what we view fe.g.. lKorista et all 
ll997l ; lLawrencell201ll) a nd some tentative ev i dence may ex- 
ist to support this (e.g.. iRichards et al. l l201ll ; lKrucz"ek~ et al] 
2011). If this is the case, this approach is still at least ener- 
getically viable. Whether the energy holds the same shape 
as when it was originally emitted from the accretion disc 
or not, it has still been emitted and should still be counted 
towards the bolometric luminosity. We note one caveat to 
this: if there is a large peak in emission that is unobserved 
in the EUV it will go unaccounted for and our calculation of 
bolometric luminosity will underestimate the true value as 
a result. If emission is reprocessed in the same wavelength 
regime, double counting will not be an issue. Uncertainty 
will still come from how isotropically different physical re- 
gions emit. 

Integrated luminosities assuming isotropy (Li ao ) are 



presented in Table [T 
ble, is calculated by 
GHz)/f(4215 A). 



Radio-loud n ess, g iven in that ta- 
Shang et al] (|201ll ) from R* = f(5 



4 DERIVING BOLOMETRIC CORRECTIONS 

Given measurements of a monochromatic luminosity and 
our estimates of integrated luminosity, it is possible to ex- 
plore a variety of bolometric corrections. We made mea- 
surements of monochromatic luminosity at three commonly 
used wavelengths, 1450, 3000, and 5100 A. These wave- 
lengths correspond to the optical window for low, medium, 
and high-redshift objects, respectively. When deriving opti- 
cal/UV corrections we did not distinguish between RL and 
RQ objects, which we justify in section BUI For bolometric 
corrections to the 2-10 keV luminosity see section [5] 



4.1 Linear Bolometric Corrections 

Historically, bolometric corrections assume a linear relation- 
ship between bolometric and a monochromatic luminosity: 



Liso — C , 



(4) 
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Table 1. Physical Data 



Object Redshift log(R*)^ QA.opt log (£; so )_ log (L bol ) ^ Csioo C3000 C1450 



3C 215 


0.4108 


3.37 


-0.69 


45.77 


45.64 


6.77 


6.66 


5.54 


3C 232 


0.5297 


2.87 


-1.55 


46.36 


46.23 


4.30 


2.64 


3.06 


3C 254 


0.7363 


3.71 


-1.41 


46.17 


46.04 


6.37 


4.26 


3.47 


3C 263 


0.6464 


3.00 


-1.96 


46.81 


46.68 


6.33 


4.00 


3.00 


3C 277.1 


0.3199 


3.52 


-1.51 


45.61 


45.48 


8.74 


5.73 


3.69 


3C 281 


0.6017 


3.23 


-1.00 


46.23 


46.10 


8.04 


5.05 


3.54 


3C 334 


0.5553 


3.11 


-1.75 


46.44 


46.32 


7.30 


3.79 


2.79 


3C 37 


0.6661 


3.74 


-3.65 


46.18 


46.05 


19.54 


8.82 


7.86 


3C 446 


1.4040 


4.34 




47.01 


46.89 


3.81 


4.92 


6.58 


3C 47 


0.4250 


3.82 


-1.22 


45.97 


45.85 


12.31 


8.27 


5.05 


4C 01.04 


0.2634 


3.41 


0.15 


45.44 


45.31 


5.60 


6.39 


10.1 


4C 06.69 


1.0002 


3.32 


-0.97 


47.30 


47.18 


7.44 


5.84 


3.71 


4C 10.06 


0.4075 


2.67 


-1.74 


46.48 


46.36 


10.26 


6.01 


4.34 


4C 11.69 


1.0370 


3.78 




47.04 


46.91 


4.53 


4.18 


3.37 


4C 19.44 


0.7192 


3.42 


-0.53 


46.91 


46.79 


8.30 


5.67 


4.11 


4C 20.24 


1.1135 


3.62 


-1.41 


46.92 


46.79 


9.23 


6.37 


5.04 


4C 22.26 


0.9760 


3.26 


0.42 


46.59 


46.47 


11.94 


8.68 


5.86 


4C 31.63 


0.2952 


2.93 


-2.18 


46.61 


46.49 


10.14 


4.17 


2.54 


4C 34.47 


0.2055 


2.62 


-1.94 


46.08 


45.96 


12.92 


5.98 


6.47 


4C 39.25 


0.6946 


3.65 


-1.74 


46.91 


46.79 


8.76 


6.14 


4.54 


4C 40.24 


1.2520 


3.94 




46.60 


46.47 


7.24 


5.82 


4.32 


4C 41.21 


0.6124 


2.91 


-1.77 


46.75 


46.63 


13.01 


6.35 


3.03 


4C 49.22 


0.3333 


3.35 


-1.00 


45.99 


45.86 


14.12 


8.07 


6.04 


4C 55.17 


0.8990 


3.74 


-0.45 


46.46 


46.33 


5.54 


4.95 


4.21 


4C 58.29 


1.3740 


2.65 




47.19 


47.07 


5.98 


4.56 


2.89 


4C73.18 


0.3027 


3.20 


-0.77 


46.40 


46.27 


7.09 


3.86 


3.95 


B2 0742 


0.4616 


2.77 


-1.73 


46.46 


46.34 


4.91 


3.23 


3.76 


IRAS F07546+3928 


0.0953 


-0.6 


-1.38 


45.43 


45.31 


3.75 


2.79 


4.72 


MRK 509 


0.0345 


-0.2 


-1.64 


45.36 


45.23 


14.74 


7.86 


6.97 


OS 562 


0.7506 


3.35 


-0.81 


46.74 


46.62 


8.31 


5.80 


3.41 


PG 0052+251 


0.1544 


-0.4 


-1.73 


45.94 


45.81 


11.36 


6.35 


4.60 


PG 0844+349 


0.0643 


-1.1 


-1.02 


45.31 


45.18 


8.35 


4.91 


4.50 


PG 0947+396 


0.2057 


-0.5 


-1.29 


45.79 


45.67 


10.14 


5.42 


4.42 


PG 0953+414 


0.2338 


-0.2 


-2.16 


46.42 


46.29 


10.80 


5.23 


3.50 


PG 1001+054 


0.1603 


0.05 


-2.12 


45.13 


45.01 


5.37 


3.30 


2.98 


PG 1100+772 


0.3114 


2.65 


-2.04 


46.46 


46.33 


10.45 


6.16 


3.64 


PG 1114+445 


0.1440 


-0.9 


-1.50 


45.39 


45.26 


3.91 


2.82 


4.34 


PG 1115+407 


0.1541 


-0.4 


-1.36 


45.71 


45.59 


12.82 


6.69 


4.06 


PG 1116+215 


0.1759 


-0.1 


-2.11 


46.31 


46.19 


16.42 


7.86 


4.25 


PG 1202+281 


0.1651 


0.04 


-2.77 


45.41 


45.29 


8.22 


5.19 


8.59 


PG 1216+069 


0.3319 


0.67 


-1.91 


46.31 


46.18 


6.66 


6.44 


4.89 


PG 1226+023 


0.1576 


3.22 


-2.00 


46.96 


46.83 


9.88 


5.05 


3.23 


PG 1309+355 


0.1823 


1.38 


-1.38 


45.74 


45.62 


4.73 


3.55 


3.42 


PG 1322+659 


0.1684 


-0.7 


-1.86 


45.74 


45.61 


11.06 


6.66 


4.64 


PG 1351+640 


0.0882 


0.09 


-1.19 


45.30 


45.18 


2.55 


3.31 


3.99 


PG 1352+183 


0.1510 


-0.6 


-1.84 


45.61 


45.48 


12.32 


6.40 


4.44 


PG 1402+261 


0.1650 


-0.5 


-2.14 


46.02 


45.90 


15.28 


6.54 


4.12 


PG 1411+442 


0.0895 


-0.8 


-1.80 


45.33 


45.20 


4.57 


2.77 


2.96 


PG 1415+451 


0.1143 


-0.5 


-1.08 


45.37 


45.24 


8.25 


6.03 


4.63 


PG 1425+267 


0.3637 


2.31 


-1.84 


46.05 


45.92 


5.88 


3.33 


3.21 


PG 1427+480 


0.2203 


-1.5 


-1.82 


45.79 


45.66 


11.24 


5.91 


3.78 


PG 1440+356 


0.0773 


-0.7 


-1.44 


45.64 


45.51 


10.91 


6.45 


4.68 


PG 1444+407 


0.2673 


-1.0 


-2.17 


46.19 


46.07 


10.71 


4.85 


4.41 


PG 1512+370 


0.3700 


2.85 


-1.97 


46.22 


46.09 


11.57 


4.88 


4.62 


PG 1543+489 


0.4000 


0.13 


-2.15 


46.30 


46.17 


8.33 


4.49 


3.15 


PG 1545+210 


0.2642 


3.00 


-1.69 


45.93 


45.80 


7.178 


4.39 


3.46 


PG 1626+554 


0.1317 


-1.0 


-1.92 


45.68 


45.55 


11.70 


5.94 


4.15 


PG 1704+608 


0.3730 


2.82 


-1.05 


46.49 


46.36 


5.89 


4.78 


3.82 


PG 2349-014 


0.1740 


2.74 


-1.48 


45.77 


45.65 


9.98 


5.56 


4.80 


PKS 0112-017 


1.3743 


3.45 




46.88 


46.76 


6.46 


4.41 


2.80 


PKS 0403-13 


0.5700 


3.73 


-0.94 


46.36 


46.24 


6.78 


5.24 


4.67 


PKS 1127-14 


1.1870 


3.88 


-1.88 


47.15 


47.02 


8.70 


6.29 


4.50 


PKS 1656+053 


0.8890 


3.10 


-1.09 


47.05 


46.93 


6.29 


5.51 


4.12 



a Radio loudness, R = f(5 GHz)/f(4215 A) 
° Integrated luminosity assuming isotropy. 

c Integrated luminosity adjusted for viewing angle: Lj, ; w 0.75 L{ ao 
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Figure 3. Bolometric luminosities resulting from y 1 / 3 , .Matthews & Fcrland il987T ). and iKorista et aT. I J1997I1 interpolation methods 
plotted against bolometric luminosities resulting from the direct power law interpolation that is ultimately employed. Solid red lines are 
best-fitting models to the data, black dashed lines are one-to-one lines. 



where L\ is the monochromatic luminosity and the factor £ 
is the bolometric correction. Using a bolometric correction 
of this form assumes a single SED shape scaled to match 
a monochromatic luminosity. Quasar SEDs vary in shape 
and there is no reason to assume a direct linear relationship 
with a zero intercept, although available data has rarely sup- 
ported more detailed fittings. 

The ratios of Li SO to A La for individual objects are 
given in Table [1] Fig. [4] shows the distribution of values 
for this ratio in each wavelength window. The 'FUSE-HST' 
subsample is known to be UV-bright and may introduce a 
slight bias towards higher UV luminosities. This would have 
the effect of reducing dispersion around the mean at 1450 
A, although it is not solely responsible for this behavior. 

We found a similar behavior of bolometric correction 
distrib utions with wavelength to iNemmen fc Brothertonl 
(2010). Bolometric correction distributions become tighter 
and shift to a lower mean with decreasing wavelength (see 
Fig. [4} . This is due to variations in SED shape from object 
to object. The peak of the Big Blue Bump is closer 1450 
A which makes bolometric correctio ns there less sensitive to 
variat ions in the shape of the bump. lNemmen fc Brothertonl 
(2010) find mean bolometric corrections (included in Table 
[2]) of 3, 5.9, and 7.6 for 1450, 3000, and 5100 A respectively. 
Note that thei r limits of integration are IR (3 /im) to X-ray 
(414 eV). The lHubenv et alJl^OOOl ) models drop off in the 
IR and are orders of magnitude be low the peak so a large 
differe nce does not result from the lNemmen fc Brothertonl 
|2010l ) lower limit compared to ours. 

Finally, we provide best-fitting bolometric corrections 
for Li 30 versus A L\ . These and other fits in this paper are 
made by minimiz ing the chi-squared statistic using mpfit 
(|Markwardtll2009l ) which employs the Levenberg-Marquardt 
least squares method. The uncertainties in A L\ are con- 
sidered to be negligible compared to those in Li ao because 
the spectra are high signal-to-noise and the monochromatic 
fluxes are well known. There are many sources of uncertainty 
in Li 30 including, but not limited to, the interpolation, varia- 
tion in SED shape, host galaxy corrections, and orientation. 
While there is uncertainty in the integration process asso- 
ciated with our selection of a model in NIR and X-ray as 



discussed, we avoid more robust fitting methods and choose 
only to minimize the chi-squared statistic. 

The resulting linear bolometric corrections are 4.2±0.1, 
5.2±0.2, and 8.1±0.4 at 1450, 3000, and 5100 A, respec- 
tively. The uncertainties given here are the 1-a errors in the 
fit coefficients. These are listed in Table [2] and compared to 
other values from the literature. 

There are several differences between the bolometric 
corrections derived here and t h ose derived in previou s stud- 
ies (namely lElvis et all fl99i I Richards et al.ll2006l ). Both 
integrated the full SED, Richards et al. ( 20061 ) specifically 
from 100 /im to 10 keV, in order to determine the bolometric 
luminosity. Here we integrate o n ly the NIR -through- X-ray 
region . As a result, lElvis et aU (1 19941 ) and iRichards et al.l 
(2006) typically find larger bolometric corrections. 

Figures [5] [6] and [7] illustrate the difference between 
inte gration methods e mplo yed here and in the literature. 
The lElvis et~ai1 l| 19941 ) and [Richards et all (|2006h bolomet- 
ric corrections will systematically overestimate bolometric 
luminosity compared to our correcti o ns du e to their inclu- 
sion of IR emission. The lElvis et alj (| 19941 ) linear relations 
give bolometric luminosities that are a factor of 1.22 1.19, 
and 1.54 greater than th ose derived here at 14 50, 3000, and 
5100 A, respectively. For IRichards et all (|2006l ) we find that 
their linear relations return bolometric luminosities a factor 
of 1.08 and 1.27 larger than those derived here at 3000 and 
5100 A, respectively. 

In ord e r to make a self-consistent com parison to the 
lElvis et al.l (| 19941 ) and I Richards et all ([2006) corrections, we 
obtained their data and reintegrated their composite SEDs 
using our limit s of integration. S ince we only integra t ed the 
composites of lElvis et alj (1 19941 ) and IRichards et all (|2006l ) 
and did not determine bolometric corrections for individ- 
ual objects, we cannot compare directly to our bolometric 
corrections from Table [2] Instead we derived bo lometric cor- 
rectio ns from the RQ composite spectrum from lShang et all 
to make our co mpari son. When we d i d this , we found 
the lElvis et al.1 (| 19941 ) and IRichards et al.l |2006l) relations 
will rather underestimate bolomet ric luminosity. Co rrections 
derived from the RQ composite of lElvis et al.l (|l994T ) will un- 
derestimate bolometric luminosity compared to corrections 
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Figure 4. Histograms showing distribution of individual linear bolometric corrections with zero intercepts at 1450, 3000, and 5100 



Table 2. Bolometric Corrections for (Li so ) = CAL;, 



Reference 


Cl450 


C3000 


C5100 


Number 
of sources 


Range in 
log(Li, oi ) 


Standard error in mean 
for 1450/3000/5100 A 


Elvis et al. (19941 a 


5.12 


6.19 


12.45 


47 


44.86-46.92 




Recalculated Elvis et al. ( 1994) b 


3.15 


3.82 


7.68 








Richards et al. (2006) 




5.62 


10.33 


259 


45.06-47.43 


••■ /0.07/0.13 


Recalculated Richards et al. (2006) c 


2.33 


3.11 


5.53 








Nemmen & Brotherton (2010) d 


3.0 


5.9 


7.6 


280 


44.60-48.50 


0.3/0.8/1.9 


This Work 


4.2 


5.2 


8.1 


63 


45.13-47.30 


0.1/0.2/0.4 



a Recalculated from the full lElvis et al] l|l994l ) radio-quiet mean SED for 1450, 3000, and 5100 A. 
^ Corrections are for Li so recalculated for the I El vis et al ] jl994t ) radio-quiet mean SED from 1 mil 
to 8 keV. 

c Corrections are for L; so recalculated for the lRichards et al. (2006) mean SED containing all objects 
from 1 fim to 8 keV. 

^ These corrections are based on lHubenv et aL I ll2000l ) theoretical accretion disc models which are 
known to produce insufficient X-rays that drop off steeply before 0.1 keV. 



derived from the RQ composite o f lShang et all (|201ll ) at all 
luminosities by factors of 0.75, 0.76, and 0.90 at 1450, 3000, 
and 5100 A, respectively . Corr ections derived from the RQ 
composite of lElvis et al] |l994h will underestimate bolomet- 
ric luminosit y compared t o corre ctions derived from the RQ 
composite of IShang et al] (|201 lh at all luminosities by fac- 
tors of 0.55, 0.62, and 0.65 at 1450, 3000, and 5100 A, respec- 
tively. The differences in these bolometric corrections reflect 
the differences in the SEDs from which they are derived and 
highlight the importance of sample s election in d e rivin g and 
applying bol ometric corrections. The lElvis et al] (|l994 ) and 
similarly the lShang et al.l l|201ll ) composites are from sam- 
ples that are known to be UV br ight and thus require l arger 
bolometric correct ions than the lRichards et all (|2006t ) . See 
IShang et alj (|201ll ) figures 6 and 8 for a comparison of the 
composites. Generally, it appears that the difference in bolo- 
metric corrections is the result of different samples as well 
as a different method and integration limits. 



4.2 Linear Bolometric Corrections with Non-Zero 
Intercepts 

W e derived a second bolometric c orrection of the form used 
bv lNemmen fc Brotherton! l|2010l ): 



log(L ls 



A + Blog(AZo, 



(5) 



We fit a linear relationship to the luminosity data in 
log-space with both A and B allowed to vary freely. We 
found that, in general, this improved the reduced \ 2 from 
the linear fit with a zero intercept. 

At 1450 A, adding a nonzero intercept significantly 
improved the fit. Using the regression analysis package in 
MiniTab® Statistical Software q we found that the proba- 
bility that the best-fitting line has a zero intercept is 4x 10~ r 
and the t-ratio is 4.71. At 3000 A the fit is only marginally 
improved by allowing the intercept to vary from zero. The 
probability that the best-fitting line has a zero intercept is 
0.151 and the t-ratio is 1.45. At 5100 A adding a nonzero 
intercept to the fit is an improvement. The probability that 



.minitab.com 
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the best-fitting line has a zero intercept is 0.004 and the 
t-ratio is 2.95. 

In practice, this provides a moderate improvement over 
using other bolometric corrections. There is no systematic 
over or underestimate for all luminosities between these 
and our bolometric corrections with a zero intercept. Ex- 
cept for the 145 A correcti o n at the lowest luminosities, 
corrections from lElvis et al.l (1 19941 ) that include IR emis- 
sion will overestimate bolometric luminosity compared to 
our bolometric corrections with nonzero intercepts by fac- 
tors of 1.08, 1.16, and 1.45, for log(AL A ) = 45.0 at 1450, 
30 00, and 5100 A , respe ctively. The bolometric corrections 
of iRichards et all (|2006T > that include IR emission will also 
overestimate bolometric luminosity compared to our bolo- 
metric corrections with nonzero intercepts by factors of 1.05 
and 1.20 for log(A£ A ) = 45.0 at 3000 and 5100 A, respec- 
tively. 

The approach th a t we employ is like that of 
iNemmen fc Brothertonl (|2010h whose fits are tabulated 
alongside ours in Table [3] They choose similar limits of in- 
tegration for calculating bolometric luminosity and find a 
non zero intercept to be a good fit as we do here. Note that 
the lHubenv et al l (2000) accretion disc models are theoret- 
ical and do not produc e X-rays. 

In comparison, theE emmcn fc Brothertonl (|2010h bolo- 
metric corrections show dif ferent behaviors depend i ng on 
wavelength. At 1450 A the INemmen fc Brothertonl (|2010l ) 
correction systematically overestimates bolometric luminos- 
ity compared to ours, by a factor of 1.50 at log(A La) = 45.0. 
At 3000 A their correction underestimates ours at most 
monochromatic luminosities in our range, by a factor of 
0.92 at log(ALx) = 45.0. At 5100 A their correction under- 
estimates ours at most monochromatic luminosities in our 
range, by a factor of 0.93 at log(AL A ) = 45.0. 

Figs. [5] [6] and show the relationship between Li SO 
and \L\ in log-log space for A = 1450, 3000, and 5100 A 
respectively. Bolometric corrections fro m the literat u re are 
show n for comparis o n. Co rrections from El vis et all (jl994h 
and IRichards et al. (2006) are those derived from the full 
SED including the IR, although those of lElvis et all (119941 ) 
have been corrected to our fiducial wavelengths. 

We calculated isotropic and bolometric luminosity for 
the 22 objects with incomplete wavelength coverage using 
log(L iso ) = (4.74±1.00) + (0.91±0.02) log(1450 L1450) from 
Tabled These values are presented along with other quan- 
tities of interest in Table [4] 

4.3 First-Order Bolometric Corrections 

The distributions of bolometric corrections in Fig. [4] and 
the SEDs themselves motivate an additional term in the 
bolometric correction. The distribution at 1450 A is much 
narrower than the distribution at 5100 A. The shape and 
wavelength of the bump vary less from object to object near 
the peak of the bump than on the sides so bolometric cor- 
rections at 1450 A, which are closer to the peak, see less 
variation than bolometric corrections at 5100 A. 

In order to decrease the scatter in the 5100 A relation- 
ship, we included optical slope in a first order bolometric 
correction. Optical slope helps to account for differences 
in SED shape by describing the relationship between the 
long- wavelength side and peak of the Big Blue Bump. Op- 




log(XL Jk ergs S" 1 ) 

Figure 5. Relation between Li ao and \L\ in log-log space at 
1450 A. The solid red line shows the \ 2 fit with a nonzero in- 
tercept, the dashed-dotted black line shows our linear best-fitting 
bolometric correction with a zero intercept, the dashed-dotted 
gray line shows the linear re lation with a nonzero intercept of 
INemmen &c Brothertonl ll2010h (NB10), and the dashed g ray line 
shows the linear bolometric correction o f Elvis et al. '(1994) (E94). 
RL objects are plotted as dark blue squares and RQ objects are 
plotted as open black circles. 

tical slope is not independent of other parameters (e.g. ori- 
entation and Edding ton fraction) in accretion disc models 
l|Hubenv et al.l [20001 so the dependence of bolometric cor- 
rections on optical slope may in fact reflect these other prop- 
erties to some extent. 

We measured optical slope in most cases as a power law 
continuum at two points: 3981 A and 6310 A. In all cases, 
the measurement was taken under the H/3 line and avoided 
the region contaminated by Fe II emission. We defined <Xx,opt 
such that F\ oc A QA opt . Because of limited wavelength cov- 
erage, optical slope could only be measured in 58 of the 63 
objects with bolometric luminosities. Fig. [8] shows the dis- 
tribution of optical slopes for the full sample and only the 
RL portion. The distributions of optical slope for the full 
sample and RL and RQ subsamples are statistically similar. 
A K-S test on the full and RL sample gives Dks = 0.12 and 
a probability of 0.86 that the two distributions come from 
the same parent sample. A K-S test on the full and RQ dis- 
tributions give Dks = 0.19 and a probability of 0.54 that 
the two distributions come from the same parent sample. A 
K-S test on the RL and RQ samples gives Dks = 0.31 and 
a probability of 0.10 that the two distributions come from 
the same parent sample. 

Optical slope is provided for objects with incomplete 
SEDs in Table [4] These are not used to find the first order 
bolometric correction and are not included in Fig. [5] 

Using MiniTab® Statistical Software, we performed a 
multiple regression analysis with 5100 L5100 and a A , op t and 
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Table 3. Bolometric Corrections for log(Li so ) = A + B log(AL A ) and First-Order Bolometric Correction 



Wavelength 


Bolometric Correction (This Work) a 


Ncmmen & 


; Brotherton (2010) Bolometric Correction 










1450 A 
3000 A 
5100 A 


log(L iso ) = (4.74 ± 1.00) + (0.91 ± 0.02) log(AL A ) 
log(L iso ) = (1.85 ± 1.27) + (0.97 ± 0.03) log(AL A ) 
log(L iso ) = (4.89 ± 1.66) + (0.91 ± 0.04) log(AL A ) 


log(L iso ) 
log(L iso ) 
log(L; so ) 


= (6.70 ±0.69) - 
= (9.24 ±0.77) - 
= (11.7 ±0.93) - 


h (0.87 ± 0.020) log(AL A ) 
h (0.81 ±0.020) log(AL A ) 
^(0.76 ±0.020) log(AL A ) 


5100 A 


log(L, so ) = (1.02 ± 0.001) log(5100 L 5 i 00 ) - (0.09 ± 0.03) a A , opt 









a Uncertainties for all fits in this table are 1 <r errors in the fit coefficients. 

Table 4. Bolometric Luminosities for Other Shane ct al . (2011) Objects 



Object Redshift log(R*) ^ a x , opt log(L Iao )^ log(t u )l 



3C 110 


0.7749 


2.68 


-2.57 


46.99 


46.86 


3C 175 


0.7693 


2.99 


-1.46 


46.91 


46.78 


3C 186 


1.0630 


3.33 


-2.19 


46.66 


46.54 


3C 207 


0.6797 


3.63 


-1.45 


46.28 


46.16 


3C 288.1 


0.9631 


3.42 


-1.87 


46.61 


46.48 


4C 12.40 


0.6836 


3.03 


-0.67 


46.17 


46.04 


4C 30.25 


1.0610 


2.92 




46.27 


46.14 


4C 64.15 


1.3000 


3.37 




46.73 


46.60 


B2 1351+51 


1.3260 


2.95 




46.66 


46.54 


B2 1555+33 


0.9420 


2.99 


-2.09 


46.21 


46.09 


B2 1611+34 


1.3945 


3.76 


-1.58 


47.06 


46.93 


MC2 0042+101 


0.5870 


2.92 


-0.63 


45.68 


45.56 


MC2 1146+111 


0.8614 


2.55 


-1.01 


46.29 


46.16 


MRK 506 


0.0428 


0.49 


-1.37 


44.67 


44.54 


PC 1103-006 


0.4234 


2.94 


-1.63 


46.30 


46.17 


PC 1259+593 


0.4769 


-1.70 


-1.80 


46.83 


46.71 


PC 1534+590 


0.0303 


0.14 


-0.69 


44.43 


44.30 


PC 2214+139 


0.0657 


-1.40 


-1.51 


45.55 


45.42 


PC 2251+113 


0.3253 


2.46 


-1.62 


46.35 


46.23 


PKS 0859-14 


1.3320 


3.43 


-1.29 


47.22 


47.10 


PKS 2216-03 


0.8993 


3.23 


-0.47 


46.95 


46.82 


TEX 1156+213 


0.3480 


2.38 


-1.51 


46.01 


45.89 



a Radio loudness, R = f(5 GHz)/f(4215 A) 

^ Calculated luminosities from log(Li so ) = (4.74 ± 1.00) - 

c Calculated luminosity adjusted for viewing angle: L\y i 



(0.91 ±0.02) log(1450Li4 50 ). 

2 0.75 Li so 



found that adding optical slope is significant with a t-ratio 
of 2.55. The probability that optical slope coefficient is zero 
is 0.014. The bolometric correction is: 

log(L lso ) = (1.02+0.001) log(5100L 5 ioo)-(0.09±0.03) a x>opt .(6) 

The uncertainties given here are standard deviations in 
the coefficients. This bolometric correction, also listed in Ta- 
ble|31 gives bolometric luminosities larger than those derived 
fr om other 5100 A c orrections given here, sm aller than those 
of Elvis et ail |l994;) , and similar to those of iRichards et all 
(2006). 

There are two potential concerns with including optical 
slope in a bolo metric correction: a break in optical slope, 
and reddening. IVanden Berk et alj (|200ll ) find that optical 
slope breaks near 5000 A in composites of SDSS quasars and 
conclude that this is partially due to contamination from 
the host galaxy and partially a real change in the contin- 
uum. If this is present in our objects, it might cause us 
to overestimate Li SO by including some light from the host 



galaxy. IVanden Berk et al.l l|200ll) find an optical slope of 
a A = —1.56 for their geometric mean composite spectrum 
over the wavelength range of 1300-5000 A (shortward of the 
break), and a slope of a A = 0.45 after the break. We found 
a mean value for optical slope of a Aj0p t = —1.53 ± 0.65 so it 
seems likely that this issue is minor in our sample. 



Bolometric corrections based on optical slope are sub- 
ject to the effects of intrinsic reddening in the SEDs. If 
present, reddening would cause more positive optical slopes, 
which would in turn underestimate true bolometric lumi- 
nosity. We have visually inspected the SEDs and do not see 
evidence of intrinsic reddening except perhaps in a few ob- 
jects. We do not think that reddening is significant in this 
sample, but it has not been systematically removed from the 
SEDs which should be noted when using these corrections. 
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Figure 6. Relation between L; so and \L\ in log-log space at 3000 
A. The solid red line shows the \ 2 fit with a nonzero intercept, the 
dashed-dotted black line shows our best-fitting bolometric correc- 
tion with a zero intercept, the dashed-dot ted gray line shows the 
linear relation with a nonzero intercept of lNemmen fe Brothertonl 
I^Oldf ) (N B10), the dotted gray line shows the bolometric correc- 
tion from iRichards et alj {2006) (R06), and the dashed g ray line 
shows the linear bolometric correction of lElvis et al.l Jl994l ) (E94). 
RL objects are plotted as dark blue squares and RQ objects are 
plotted as open black circles. 



Figure 7. Relation between Li ao and \L\ in log-log space at 
5100 A. The solid red line shows the \ 2 fit with a nonzero in- 
tercept, the dashed-dotted black line shows our linear best-fitting 
bolometric correction with a zero intercept, the dashed-dotted 
gray line shows the linear re lation with a nonzero intercept of 
iNemmen fc Brothertonl J2010f) (NB10), the dot t ed gra y line shows 
the bolometric correction from [Richards et al, (2006) (R06), and 
the dashed gray line shows the linear bolometric correction of 
lElvis et al] jl994h (E94). RL objects are plotted as dark blue 
squares and RQ objects are plotted as open black circles. 



4.4 Isotropy and Viewing Angles 

Based on studies of theoretical accretion disc models, we 
suggest a correction to bolometric luminosity that results 
from our assumption of isotropy. An accretion disc viewed 
face-on (i.e. at 0°) will be UV-brighter than one viewed 
edge-on. Bolometric corrections based on face-on observa- 
tions may thus overestimate bolometric luminosity. There 
are intermediate angles where isotropy is a reasonable as- 
sumption (theoretically near 55° according to fig. 11 of 
INemmen fc Brothertonl [2010), but quasars may not be 
viewed at those angles on average. 

There is evidence to suggest that, at a viewing angle 
of 55°, we would not be able to see into the central region 
of the AGN to view a quasar and a c o rrecti on for assuming 
isotropy is wa rranted. iNenkova et all {2008) compare dusty 
torus models {Nenkova et al. 20081 ) to observations of type 
1 and type 2 sources and find general agreement for opening 
angles in the range 40° — 60°. For torus opening angles in 
this range, a line of sight at 55° would intersect the dusty 
torus. 

If the viewing angle is known, or can at least be esti- 
mated, it would be possible to make a correction for the 
assum ption for isotropy based on viewing angle. iBarthell 
{l989( ) estimates an average viewing angle from his study o f 
a sample of 3CRR {Laing et al.lll983l ; ISpinrad et al.lll985l ) 
bright radio sources with 0.5 < z < 1.0 selected at 178 



MHz. At low radio frequencies, emission is dominated by 
optically thin synchrotron radiation from radio lobes which 
emit isotropically, leading Barthel to assume that the sam- 
ple was not significantly biased with respect to orientation. 
He used the ratio of broad-lined quasars to narrow-lined ra- 
dio galaxies, presumably quasars obscured by a torus along 
our line of sight, to estimate an average viewing angle of 
quasars to be 31°. This angle is determined by integrating 

r-44 4° 

J 0O ' 9sin(0) d9, where the limits of integration are a face- 
on disc and the edge of the dusty torus. 

If we assume that the IShang et al] {201lT ) sample is 
randomly selected and is preferent i ally v iewed at the an- 
gle of 31° suggested by IBarthell l|l989l ). bolometric lu- 
minosity will be over-estim ated by about 33 per cent 
(N cnimcn fc Brothertonll2010l . fig. 11.). We therefore suggest 
that the bolometric corrections derived here be decreased as 
follows: 



Lboi — J Li, 



' 0.75 Liso 



(7) 



where / is a factor describing the bias of an average view- 
ing angle due to the anisotropic emission from a disc. For 
objects to which the line of sight is known, a theoretical 
correction as a function of vie wing angle can be found from 
INemmen fc Brothertonl (|2010l ). fig. 11 and applied to Equa- 
tion [6] Their predicted range for this correction to Li SO is 
approximately 67 to 200 per cent for a relativistic disc in- 
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Figure 8. Histogram of optical slopes measured in 58 objects 
with complete wavelength coverage. The shaded histogram shows 
the distribution of optical slopes for 35 RL objects. A K-S test 
gives Dxs = 0.12 and a probability of 0.86 that the two distri- 
butions are from the same parent sample. 



eluding projection effects and limb darkening at viewing an- 
gles of zero to ninety degrees. 

We remind the reader that our sample may not repre- 
sent a random distribution of viewing angles. The 'FUSE' 
and 'PGX' subsamples were selected to be UV-bright. This 
means those subsamples may be biased to more face-on ob- 
jects (i.e. smaller viewing angles). The 'RLQ' sample may be 
biased towards more edge-on orientations (i.e. larger view- 
ing angles) because the blazars (which are viewed pole-on) 
were removed. 

Since the relations of lElvis et all l]l994h and 
iRichards et al.l (|2006l ) already overestimate bolometric 
luminosity compared to corrections provided here, applying 
the correction for viewing angle further widens the gap. 
Adopting the correction for isotropy in addition to one of 
our bolometric corrections will give bolometric luminosi- 
ties as much as a f actor of two lowe r, on average, than 
those resulting from lElvis et alj fl994) and [Richards et al.l 
(2006) corrections. As before, the difference in bolometric 
luminosity at 1450 and 3000 A is less than at 5100 A, and 
bolometric luminosit ies resulting from the corrections of 
IRichards etafl (2006) are more similar to ours than those 
resulting from the corrections of [Elvis et al l l| 19941 ). 



4.5 Higher Order and Emission-Line Bolometric 
Corrections 

We might expect that a higher order bolometric correction 
may fit the data due to the relationshi p between optic al- X- 
ray slope, a ox , and luminosity (e.e.. I Just et al.lfeoO'it ). We 
investigate the data and see no need for a higher order bolo- 



metric correction over this range at these wavelengths. It 
is possible that at higher luminosities the data will depart 
from a linear model. 

It may be possible to provide more accurate bolo- 
metric corrections based on correlations between emission- 
line properties and properties of the SED. For example, 
iKruczek et"afl (|201lT ) measure C IV equivalent width (EW) 
and velocity shift (relative to the quasar rest frame) for a 
sample of radio-quiet SDSS quasars and determine that ob- 
jects that fall on opposite ends of the EW-velocity shift con- 
tinuum will require different bolometric corrections. Specif- 
ically, bolometric corrections will be underestimated for ob- 
jects with large EW and small blueshift and overestimated 
for objects with small EW and large blueshift. Emission-line 
bolometric corrections of this type are beyond the scope of 
this paper and may be pursued in the future. 

4.6 Radio Class Issues 

In the previous sections we have not differentiated between 
RL and RQ objects when deriving optical/UV bolometric 
corrections. This treatment requires specific justification be- 
cause RL and RQ quasars are known to have differences in 
their SEDs. 

T he behavior of the SEP at X-ray wavelengths is well 
known. fZamorani et al.l l|l98 ll ) find that RL quasars are typi- 
cally three times more X-ray lum inous than RQ quas ars with 
similar optical/UV luminosities. iMiller et all (|201li ) charac- 
terize the X-ray excess of RL quasars with respect to RQ 
quasars and find that it can be as high as a factor of 3.4- 
10.7 for extremely radio-loud quasars with R* > 3.5. These 
differences might indicate that different bolometric correc- 
tions as a function of radio class may be appropriate. We 
did not assume a priori that RL and RQ quasars can be 
combined. 

We made fits of bolometric luminosity versus monochro- 
matic luminosity to determine whether RL and RQ re- 
quire separate bolometric corrections. Because the RL ob- 
jects have higher luminosities on average, we isolated radio- 
loudness by fitting only the objects in the region where RL 
and RQ overlap in luminosity. 

Fig. El where panels (a), (b), and (c) are for 1450, 3000, 
and 5100 Arespectively, shows the full sample with RL ob- 
jects in the overlapping region plotted as blue squares and 
RQ objects in that region plotted as open black circles. The 
best-fitting lines are solid blue and dashed-dotted black re- 
spectively. The dashed lines show the 95% confidence inter- 
vals which are shaded blue for RL and cross-hatched in gray 
for RQ. 

Based on Fig. [9] we determined that the bolometric cor- 
rections for RL and RQ are consistent within the 95% con- 
fidence intervals at optical/UV wavelengths. Relaxing our 
fits to include all of the objects does not change this re- 
sult, nor does comparing confidence intervals for fits to the 
RL and RQ subsamples. This is likely because, while the 
X-ray behavior in RL and RQ objects is different, it is not 
energetically significant compared to the Big Blue Bump 
(which is similar in RL and RQ quasars, see the SEDs in 
IShang et al]|2011i ). Based on this analysis, we combine RL 
and RQ quasars for optical bolometric corrections. 

The caveat to combining RL and RQ optical/UV bolo- 
metric corrections is the unknown emission in the extreme 
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UV and soft X-ray; RL and RQ quasars may require a dif- 
ferent interpolation in this region. While it is not currently 
possible to observe directly whether RL quasars have more 
ionizing flux in the EUV-soft-X-ray than RQ quasars, it may 
be possible to infer it from emission-line properties, though 
this is beyond the scope of this paper. 



5 X-RAY BOLOMETRIC CORRECTIONS 

We provide several bolometric corrections to the 2-10 keV 
X-ray luminosity, which is calculated by extrapolating the 
data to 10 keV with a power law of the same slope as the 
highest energy X-ray data and integrating the SED between 
2 and 10 keV. As with the optical/UV corrections, we first 
determined whether it is appropriate to treat RL and RQ 
objects together. 

To evaluate the 2-10 keV X-ray bolometric correction 
for RL and RQ objects we again fit them separately and 
compared the fits within the 95% confidence intervals. Fig. 
[I0]panel (a) shows the 2-10 keV bolometric correction where 
we have fit the RL and RQ objects separately in the region 
of overlapping luminosity. While the lines appear to differ, 
there are so few points in this region that the confidence 
intervals are very large. Within 95% confidence the lines 
are consistent. Due to the large uncertainty in these lines, 
we also fit the RL and RQ data without matching lumi- 
nosity range. In Fig. [10] panel (b) we see that these lines 
are not consistent within 95% confidence when luminosity is 
not matched. The difference in RL an RQ objects at X- ray 
energies has been well characterized (jMiller et al.ll201ll ) so 
when making X-ray bolometric corrections we provide sep- 
arate RL and RQ X-ray bolometric corrections. 

We made bolometric corrections of the linear form with 
and without a zero intercept fit to the entire sample, the 
RL and RQ objects separately, and the RL and RQ ob- 
jects separately in the region of overlapping luminosity. The 
fits with the nonzero intercept always had a lower chi- 
squared than those where the intercept is fixed at zero and so 
are preferred. More complex fits (e.g.. [Marconi et al.ll20oj ; 
Hop kins et ai1l2007l ) are not necessary to fit these data over 
this luminosity range. Table [5] lists X-ray bolometric correc- 
tions. 

X-ray bolometric corrections have been shown to be 
luminosity dependent ((Marconi et al.l |2004| ; iHopkins et al.l 
120071 ; IVasudevan fc Fabianl 120071 1. We investigated the de- 
pendance of our X-ray bolometric corrections on luminosity 
via a ox - This depend ence has specifically b een character- 
ized, most recently bv lMarchese et all (|201lh for 195 X-ray 
selected Type 1 AGN from the XMM - Newton serendipitous 
survey and also iLusso et al l (|201Cf ) for 545 Type 1 AGN 
from the XMM-COSMOS survey. We calcula ted a ox from 
the 25 00 A luminosity using the relation from Steffen et al.l 
|2006l ): 



a„x = -0.137 • log(L 25 oo) + 2.638. 



(8) 



Our objects have luminosities within the range for 
which this relation was derived. We found only a weak trend 
with lots of scatter of X-ray bolometric corrections with a ox , 
if any at all, and so do not include a luminosity-dependent 
bolometric correction. The lack of a strong rela t ionsh ip in 
our data where one was found bv lMarchese et all (|201ll ) and 
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Figure 11. Comp arison between our X - ray b olometric correc- 
tions and those of Vasudcvan & Fabian (2007). The fit to the 
RL data in the region of overlapping luminosity with our RQ 
objects, shown as blue squares, is given by the solid blue line 
with hashed light blue 95% confidence intervals. The fit to the 
RQ data in the region of overlapping luminosity, shown as open 
black circles, is given by the dashed-dotted black line with hashed 
gray 95% confidence intervals. Our high-luminosity RL objects 
and low-luminosity RQ ob jects are shown as gold crosses. The 
IVasudevan fc Fabianl (2007) objects are shown as green triangles. 
Ob jects of overlap betwe e n our samples have been removed from 
the IVasudevan &: Fabianl (J2007J) sample. The lone RL object that 
sits far outside of the 95% confidence intervals in the upper half 
of the plot was removed from the fit as an outlier. 



ILusso et al.l (|2010l ) may be due to weak X-ray sources that 
incr ease the scatter in our samp le. 

IVasudevan fc Fabianl (|2007l ) calculate X-ray bolometric 
corrections with a zero inte rcept for 54 AGN wit h model 
SEDs fit to optical data from Ba skin fc Laori l|2005l ) and the 
NASA/IPAC Extragalactic Database and UV/X-ray data 
from FUSE. There is some overlap between their sample and 
ours, but because quantities are calculated from different 
SEDs and limits, we found only general agreement with their 
X-ray bolometric corrections. Fig. | 11 | shows o ur data plot- 
ted with those of IVasudevan fc Fabianl (|2007l '). Their data 
and ours seem to be reverse-cor related but with large scat- 
ter. IVasudevan fc Fabianl l|2007l ) point out that the X -ray 
bolometric corrections derived bv lHopkins et al.l (2007) and 
IMarconi et al.l l|2004h seem to have inadequate dispersion to 
describe the data. Our data support this conclusion. 

We caution the reader that mean X-ray bolometric cor- 
rections are very inaccurate for individual objects due to 

2 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 
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Figure 9. Fits to bolometric luminosity versus monochromatic luminosity for RL and RQ quasars. Data are used for fitting only in the 
region where RL and RQ quasars overlap in luminosity. RL quasars are plotted as blue squares and the best-fitting line is solid blue. RQ 
quasars are plotted as open black circles. Black crosses are RL and RQ objects that do not overlap in luminosity. The best-fitting line is 
a dashed-dotted black. Dashed lines indicate 95% confidence intervals. RL quasar confidence intervals are shaded light blue, RQ quasar 
confidence intervals are shaded with gray cross-hatch. Figures (a), (b), and (c) are for 1450, 3000, and 5100 A, respectively. 



Table 5. Bolometric Corrections for 2-10 keV Luminosity 



Sample Linear w/ zero intercept Linear w/ nonzero intercept 



Full Sample 




= (38.00 


± 


06.84) L 2 - 


-lOkeV 


log(L iso ) 


= (25.14 


±01.93) - 


- (0.47 


± 


0.043) 


lot 




-lOfceVj 


Radio-Loud 


^iso 


= (23.31 


± 


03.91) L 2 . 


-lOkeV 


lot 


S(^iso) 


= (23.04 


± 03.60) - 


- (0.52 


± 


0.080) 


toe 


;(L 2 . 


-lOfceV") 


Radio- Quiet 


^iso 


= (88.99 


± 


30.18) L 2 - 


-lOkeV 


lo t 


5(^iso) 


= (33.06 


±03.17) - 


- (0.29 


± 


0.072) 


lot 


;(L 2 - 


-lOfeev) 


L-matched Radio-Loud 


L%so 


= (43.86 


± 


13.30) L 2 - 


-lOkeV 


lo t 


S(^iso) 


= (29.30 


± 18.16) - 


- (0.38 


± 


0.409) 


lot 


;(L 2 . 


-lOfceV") 


L-matched Radio-Quiet 


^iso 


= (45.56 


± 


08.38) L 2 . 


-lOkeV 


lot 


S(^iso) 


= (15.32 


± 11.69) - 


- (0.69 


± 


0.264) 


las 


;(L 2 _ 


-lOfceV") 
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Figure 10. Fits to bolometric luminosity versus X-ray luminosity for RL and RQ quasars. In panel (a), data are used for fitting only 
in the region where RL and RQ quasars overlap in luminosity. RL quasars are plotted as blue squares and the best-fitting line is solid 
blue. RQ quasars are plotted as open black circles. The best-fitting line is a dashed-dotted black. Dashed lines indicate 95% confidence 
intervals. RL quasar confidence intervals are shaded light blue, RQ quasar confidence intervals are shaded with gray cross-hatch. In panel 
(b) all of the data are used for fitting and the same conventions apply. 



variation in the X-ray emission from object to object. This 
is in agreement with studies by iMarchese et all (|201lf) who 
conclude that X-ray luminosity is not a useful proxy for bolo- 
metric luminosity due to the large intrinsic scatter in X-ray 
bolometric corrections. Fig. [12] shows the distribution of X- 
ray bolometric corrections. The bolometric corrections are 
given in log in order to display the wide range of individual 
corrections. 



6 SUMMARY 

We provide a new gene ration of bolometric corrections de- 
rived from the SEDs of IShang et~ai1 (^Ollf ). We integrated 
the NIR- X-ray region of each SED to estimate bolometric 
luminosity and derived a variety of bolometric corrections, 
although we determined that separate RL and RQ opti- 
cal/UV corrections were not warranted by the data. We also 
provided bolometric corrections to the 2-10 keV luminosity, 
including separate RL and RQ corrections. 

We provided optical/UV bolometric corrections with 
four levels of increasing sophistication and several X-ray 
bolometric corrections, listed below. These are valid under 
the accompanying assumptions. 

(i) Traditional, linear best-fitting bolometric corrections 
of the form Li SO — £ \L\ at three wavelengths. This assumes 
SED shapes do not vary from object to object and bolomet- 
ric luminosity and monochromatic luminosity are linearly 
related. 

(ii) Linear, best-fitting bolometric corrections with a 
nonzero intercept of the form log(Li SO ) = A + B\og(\L\). 

(iii) A first-order bolometric correction. The bolometric 
correction at 5100 A can be improved by adding a correction 
based on the optical continuum slope. 




12 3 4 

log(X-ray bolometric correction) 

Figure 12. Histogram of 2-10 keV X-ray bolometric corrections. 
Note that the x-axis is in log in order to display the extremely 
wide range in X-ray bolometric corrections. 



(iv) Bolometric corrections assuming isotropy should 
probably be corrected for an average viewing angle bias. 

(v) Bolometric corrections to the 2-10 keV luminosity for 
the full sample, RL and RQ subsamples, and RL and RQ 
subsamples in the region of overlapping luminosity. 
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Mean bolometric corrections derived h ere are smaller on 
average than those d erived previously fe.g. lElvis et al.|[l994l ; 
iRichards et al.l feOOrJ ) due do a difference in methodology 
and da ta. We found good agreement with the theoretical 
work of lNemmen fc Brothertonl (|2010h as the distribution of 
bolometric corrections tightens and shifts to a lower mean 
with decreasing wavelength. 

For objects with properties that fall in the region of pa- 
rameter space covered by this sample (log(L;, ;)=45.1-47.3), 
we recommend the following bolometric corrections: 

• Recommended Correction at 1450 A 

log(L lso ) = (4.74 ± 1.00) + (0.91 ± 0.02) log(1450L 145 o) (9) 

• Recommended Correction at 3000 A 

log(L lso ) = (1.85 ± 1.27) + (0.98 ± 0.03) log(3000L 30 oo) (10) 

• Recommended Corrections at 5100 A 

log(L lso ) = (4.89 ± 1.66) + (0.91 ± 0.04) log(5100L 5 ioo) (11) 
log(L lso ) = (1.02±0.001) log(5100L 5 ioo)-(0.09±0.03) a\,„ pt (12) 

• Correction for Average Viewing Angle 

This correction should be applied to bolometric luminosi- 
ties resulting from any of our other optical/UV bolometric 
corrections. 

Lboi = f Li SO ~ 0.75 Li SO (13) 

• X-ray Bolometric Corrections 

log(L lso , M ) = (23.04±03.60)+(0.52±0.080) log(L 2 _iofc e v)(14) 

log(L lso , flQ ) = (33.06±03.17)+(0.29±0.072) log(L 2 -i fc e v)(15) 

Bolometric corrections with nonzero intercepts are 
marginally preferred to those without. At 5100 A, the bolo- 
metric correction with optical slope is preferable to the one 
with the nonzero intercept. All bolometric luminosities de- 
rived here assume isotropy and we recommend they be cor- 
rected for an average viewing angle bias. X-ray bolomet- 
ric corrections are uncertain and optical/UV corrections are 
preferred when possible. 
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